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Abstract. The self-affinity of growing systems with radial symmetry, from tumors 
to grain-grain displacement, has devoted increasing interest in the last decade. In this 
work, we analyzed features about the interface scaling of these clusters through large 
scale simulations (up to 3 x 10 7 particles) of two-dimensional growth processes with 
special emphasis on the off-lattice Eden model. The central objective is to discuss 
an important pitfall associated to the evaluation of the growth exponent (3 of these 
systems. We show that the (3 value depends on the choice of the origin used to 
determine the interface width. We considered two strategies frequently used. When the 
width is evaluated in relation to the center of mass (CM) of the border, the exponent 
obtained for the Eden model was (3cm — 0.404 ± 0.013, in very good agreement with 
previous reported values. However, if the border CM is replaced by the initial seed 
position (a static origin), the exponent (3q = 0.333 ± 0.010, in complete agreement 
with the KPZ value Pkpz = 1/3, was found. The difference between (3cm and (3$ 
was explained through the border CM fluctuations that grow faster than the overall 
interface fluctuations. Indeed, we show that the exponents /?o and (3cm characterize 
large and small wavelength fluctuations of the interface, respectively. These finds were 
also observed in three distinct lattice models, in which the lattice-imposed anisotropy 
is absent. 
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1. Introduction 

Interfaces in dynamic systems are present everywhere in nature, ranging from thin film 
deposition to biological growth. The scaling analysis of these interfaces constitutes a 
procedure widely used to characterize the underlying dynamics of these growth processes 
[H EJ E] . One class of them that has attracted increasing interest along the past decade 
is the interface scaling of biological systems that exhibit radially symmetric growth 
|H El El H El El Ell EH EE] . Initially, the scaling analysis of such biological systems was 
essentially restricted to the theoretical models |U El El E] ■ However, several biological 
experiments as, for instance, plant callus evolution [8] and the growth of malignant cells 
and tumoral explants [9], [TOl [TTJ [12] , have been reported in the past few years. Self-affine 
interfaces of radially symmetric patterns were also found in experiments of grain-grain 
displacement in Hele-Shaw cells in the quasi static regime [13] . 

The Eden model [14J is the simplest discrete example that generates radially 
symmetric patterns with self-affine interfaces. It was initially designed to describe 
biological pattern formation. The original version was studied on a square lattice, 
in which an occupied site represents a cell. The simulation begins with a single cell 
at the center of the lattice and the growth rules are as follows: at each step a site of 
the cluster periphery (an occupied site with at least one empty nearest neighbor) is 
chosen at random and one of its empty nearest neighbors (NN) is selected with equal 
probability and occupied. Variants of the Eden rules were studied and the original 
model is commonly called Eden B [2j. From the biological viewpoint, the Eden model is 
unrealistic, but it produces interfaces with a rich scaling usually analyzed through the 
interface width w defined as the root mean square deviation of the interface around its 
mean value 



where a set of N distances represents the interface and f is the mean value of these 
distances. 

In the case of interfaces grown from a d dimensional substrate of linear size L pQ , 
the interface width commonly behaves as: w ~ t", for t <C L z , and w ~ L a , for t 3> L z . 
The exponents j3, a and z (growth, roughness and dynamic exponents, respectively) are 
related by a = /3z. A given set of values of these exponents defines a universality class 
and can reveal fundamental properties of the interface dynamics. Examples in 1 + 1 
dimensions with exactly known exponents include the universality classes of Edwards- 
Wilkinson (p = 1/4, a = 1/2, and z = 2) [15j, Kardar-Parisi-Zhang (J3 = 1/3, a = 1/2, 
and z = 3/2) [16], and Mullins-Herring (/3 = 3/8, a = 3/2, and z — 4) [17] . also known 
as EW, KPZ, and MH universality classes, respectively. In general, a universality class 
is related to a dominant physical process of the interface dynamics [1]. The EW, KPZ 
and MH universality classes are related to local relaxation, lateral growth, and surface 
diffusion, respectively. This interface analysis was applied to the growth of several types 
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of tumors suggests that this biological growth dynamics belongs to the MH universality 
class |10| . These experiments, which have a very significant impact because they reveal 
a universal dynamics of tumors, were grounded on a interface scaling analysis using 
Eq. (TIJ for radial tumors. However, as we will discuss along this paper, there is an 
important pitfall associated to this procedure that can lead to erroneous conclusions 
about the growth exponent (3 and, consequently, to the universality class of the process. 

Simulations of Eden clusters grown from a flat substrate show that the model 
belongs to the KPZ universality class [3j [TBI OS] • However, the shape of the Eden clusters 
grown from a seed is very sensitive to the lattice anisotropy [201 EI]- Zabolitzky and 
Stauffer [20] simulated clusters of the Eden A modejj] with N ~ 10 9 particles on square 
lattices and observed a complex behavior of the interface width. For small clusters, 
a relatively good agreement with KPZ growth exponent ((3 ~ 1/3) was observed, 
contrasting with the linear dependence on time {(3 — > 1) obtained for asymptotically 
large clusters. The value (3 = 1 is due to the diamond-like shape of the cluster imposed 
by the square lattice anisotropy [22J. In order to determine the growth exponent of 
round Eden clusters, off-lattice simulations with N « 2 x 10 5 particles were done by 
Wang et al. [23l I2T] by considering the center of mass (CM) of the cluster borders as 
the origin for the evaluation of the roughness. They measured an exponent (3 = 0.396 
and claimed that this value is close to the KPZ exponent. 

In this work, we demonstrate through large scale simulations (N > 3 x 10 7 ) that 
the border CM fluctuations in the off-lattice Eden model are not negligible and the 
random motion of the border CM determines the exponent value observed by Wang et 
al. Indeed, we found a growth exponent (3q = 0.333 ± 0.010 very close to 1/3 when an 
origin fixed on the initial seed is used. When the cluster CM (all cells are taken into 
account) is used, the exponent asymptotically converges to (3q. We also applied these 
ideas to three distinct on-lattice growth models in which the lattice-imposed anisotropy 
is absent. The paper is organized as follows. In Sec. El the simulation procedures 
for Eden model are described. In Sec. [3l the results for Eden Model are presented 
and discussed while the on-lattice models and the scaling analysis are presented in the 
section [H Finally, some conclusions are drawn in Sec. [5j 

2. Off-lattice Algorithm 

The simulations supporting the aforementioned results used the off-lattice algorithm 
proposed by Wang et al. [23], in which the particles are represented by discs of diameter 
a and the growth rules are the following. 

(i) An active cell (able to grow) is introduced on a plane. 

(ii) A cell is selected randomly among the active ones. The intervals along which an 
adjacent cell can be grown without overlapping any existing cells are identified. 

f In this version, the empty sites neighboring the cluster interface are chosen at random and then 
occupied 
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Then, a new cell is grown in a direction randomly chosen in the allowed intervals, 
(iii) If there are no possible growth directions, the cell is labeled as inactive. 

Also, since our interest is focused on the interface scaling, we introduced an optimization 
where any active cell inside a central core of radius r c is labeled as inactive. Since the 
inactivation of the particles near or belonging to the interface must be avoided, r c = 0.8f 
was chosen, where f is the mean radius of the interface. This optimization was used 
only for f > 300a. With these procedures, we grew clusters exceeding 3 x 10 7 particles 
(more than two orders of magnitude larger than those obtained in previously reported 
simulations [23]). In FIG. a typical growth pattern and the corresponding bordejj] 
are illustrated. The mean density of cells inside the clusters has the closely constant 
value p = 0.633 ± 0.001, slightly lower than the density estimated by Wang et al. 
(p = 0.6500 ± 0.0008) [23J. This difference is due to the divisions occurring inside the 
region r < r c in the original Wang algorithm which are forbidden in our modified rules. 




Figure 1. (color online) A small Eden cluster with 6000 particles. The border is 
depicted in red. 

In order to evaluate the interface width, four methods were used to define the 
distances in Eq. (CQ). In the first one, the border CM method, r« is the distance from 
the instantaneous border CM to a site on the cluster border. In the second one, the 
cluster CM method, the CM of the entire cluster is used instead of the border one. In 
the third one, the seed method, the distance from the CM is replaced by the distance 
from the initial seed (a static origin). Finally, in the last one, the sector method, the 

| The border is defined as the set of cells that forms an external layer impenetrable for incoming cells. 
Consequently, the spaces between consecutive border cells is smaller than a cell diameter. 
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border is divided in k sectors of equal angular separation and w is defined as the average 
of the standard deviation of the distances from the seed along each sector. Since the 
dynamics of the model occurs essentially on the border, we can do an analogy with the 
Eden model grown from flat substrates [3] and consider the time proportional to the 
number of peripheral particles or, equivalently, to the mean radius. 

3. Off-lattice simulations 

Figure [2] shows the interface width evaluated through the previously mentioned methods 
except the sector one. One can clearly observe distinct power laws for the roughness 
evaluated through the border CM and the seed methods. The corresponding exponents 
are @cm ~ 0.40 and f3 ~ 0.33, respectively. The first value is in very good agreement 
with the simulations performed by Wang et al. [23] (/3 = 0.396), whereas the second one 
is in excellent agreement with the KPZ universality class (Pkpz = 1/3). The last result 
confirms, for the first time, the claim that radial off-lattice Eden clusters belong to the 
KPZ universality class. The inset of FIG. EJ in which the ratios between the interface 
width evaluated through distinct procedures are plotted, shows that the cluster CM and 
seed methods have the same asymptotic scaling (f > 10 3 a or, equivalently, N > 10 7 ), 
while distinct growth exponents are observed in the intermediate intervals. Certainly, 
this transient can lead to wrong conclusions about the universality classes of experiments 
where the system can not grow forever 0, [9], EUl EH] ■ We evaluated the local slope of the 
plots In w against In f as shown in the bottom of FIG. [2j As one can see, the exponents (5q 
and Pcm are clearly distinct and oscillate around the their expected values 1/3 and 2/5, 
respectively, for f > 20a. Considering these fluctuations as estimates of the exponent 
uncertainties, we obtained j3o = 0.333 ± 0.010 and (3cm = 0.404 ± 0.013. 

The previous difference can be better understood by analyzing the CM evolution. 
Figure [3] shows three stages of typical trajectories of the cluster and border CMs along 
a simulation (top). In these walks, a step is defined as the CM displacement when the 
cluster radius of gyratiorjjj] increases by a cell diameter. The difference between the 
trajectories is evident. In the former, the CM wanders through a region with a few cell 
diameters and follows a trajectory of low fractal dimension. In the later, the CM wanders 
around a region of increasing amplitude and magnitude with the same order than the 
interface width. Consequently, the trajectory of the border CM is more compact than 
that of the cluster CM. The bottom of FIG. [3] shows the CM mean distance from the 
initial seed Rcm as a function of the mean radius. As can be seen, for both cluster 
and border CMs, the mean distance grows approximately as a power law Rcm ~ ^ 7 for 
f > 200a. The exponents 75 = 0.45 ± 0.04 and 7 C = 0.24 ± 0.05 were found for the 
border and cluster CMs, respectively. These results elucidate the differences between the 
growth exponents. When the cluster CM is used, its fluctuations around the initial seed 
grow slower than the interface fluctuations (/3 > j c ). So, in the asymptotic limit the CM 

|| The definition of radius of gyration is r g — Q^iLi rf /N) 1 / 2 , where represent the distances from 
the initial seed. 
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Figure 2. In the top, interface width evaluated using the border CM (wcm), the 
cluster CM (rop M ), and the seed (wo) methods is shown. The straight lines correspond 
to the power law fits for f > 10 2 a. The inset shows the ratios T — Wq M /wcm 
(triangles) and T = Wq M /wq (circles). In the bottom, the local slope (the local growth 
exponent) is plotted as a function of the system size. The horizontal lines represent 
the slopes 1/3 and 2/5. All these curves result from 10 3 independent samples, f and 
w are given in cell diameter unities. 

fluctuations are negligible and the growth exponent converges to (5q (inset of FIG. |2]). 
In contrast, the border CM fluctuations increase faster than the interface fluctuations 
(7b > 1/3) and, hence, they do not become asymptotically negligible. These findings 
might be easily verified in experimental essays such as those related to the callus growth 
[8], tumor evolution [91 [TOJ [11] , or grain-grain displacement in Hele-Shaw cells [13]. 

In FIG. HI the results for the sector method are confronted with those obtained 
using the border CM and seed methods. This figure shows the ratios Wk/w and Wk/wcM 
(definitions in the legend of FIG. H|) as functions of f. A very peculiar behavior arises 
from these curves: if the interface is divided in a small number of sectors, the growth 
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Figure 3. (color online) Typical trajectories of the cluster (top left) and the border 
(top right) centers of mass. The colors correspond to 342 (black), 684 (red), and 1026 
(blue) steps (definition in the text). In the bottom, the border (circles) and cluster 
(squares) CM mean distances from the initial seed are plotted as functions of the mean 
radius. Rcm and f are given in cell diameter unities. 

exponent quickly converges to whereas the scaling of border CM method is observed 
for a large number of sectors. For intermediate number of sectors (k = 12), a crossover 
from Pcm to (3q can be perceived, suggesting that the growth exponent asymptotically 
reaches the value 1/3. So, the KPZ universality class observed in the Eden clusters 
results from the large wavelength fluctuations of the interface. Also, one can infer 
that the border CM fluctuations are straightforwardly related to the small wavelength 
fluctuations of the border. This analysis may be relevant for the characterization of 
several experiments. For example, distinct scales for interface fluctuations have been 
recently identified in cell membrane of single macrophages during the phagocytosis 
process [25J, [26] 

4. On-lattice models 

In order to corroborate the scaling concepts based on the off-lattice Eden model 
simulations, we proposed three different growth rules. We chose on-lattice models due to 
the easiness for proposing and implementing new rules. However, these models should 
avoid the undesirable lattice anisotropy effects. The rules described in the sequence 
fulfill this requirement. In all models, the simulations begin with an occupied site at 
the center of the lattice and the growth rules at each time step are the following: 

(i) Model I. A particle is released at the center of the lattice (on an occupied site) and 
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Figure 4. Ratios between interface widths determined through the sector (wk), the 
seed (wo), and the border CM (wcm) methods as functions of the mean radius f. The 
corresponding numbers of sectors are indicated in the legend and r is given in cell 
diameters unities. 



follows a ballistic trajectory at a random direction while it does not reach an empty 
site (FIG.Ggieft). 

(ii) Model II. Like in Model I, except that the particle is released at any occupied site 
chosen at random (FIG. [^middle). 

(iii) Model III. Firstly, the same rule of Model II is implemented. Secondly, other particle 
is added to the symmetric position in the cluster, passing through the center of the 
lattice, with probability p (FIG. fright). This model allows one to control the CM 
fluctuations. In particular, the CM is static for p = 1. 

The isotropy of the patterns was confirmed using noise reduction methods (see [2] 
or [3] for details about the method). The growth exponents are summarized in table [TJ 
As one can see, independently of the model, the exponent (3cm has a value close to 2/5 
whereas the exponent (3q for the model I is neatly different from the other ones. This 
results together with Eden model simulations suggest that [3cm = 2/5 is a universal 
exponent. 

Model III allows one to control the CM fluctuations and, consequently, to verify 
their roles on the growth exponents. Figure [6] shows the quantities wo, wcm, Rcm and 
the ratio wcm/wq for p = 0.90 as functions of time. These curves provide evidences 
of a crossover induced by the amplification of the CM fluctuations. Indeed, the scaling 
regime with [3cm 7^ Po emerges when Rcm ~ ^o- 
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Figure 5. Schematic representation of the on-lattice models. The clusters are 
representative samples of actual simulations. Model I (left): the particle leaves the 
center towards the border; Model II (middle): the particle follows a ballistic trajectory 
with random initial positions and directions; Model III (right): a first particle is grown 
through a random ballistic trajectory (continuous line) and a second one is added to 
the opposite side of the cluster (dashed line). 



Model fa Pcm 



Model I 




0.284 


± 


0.009 


0.38 


± 


0.02 


Model II 




0.209 


± 


0.006 


0.40 


± 


0.05 


Model III (p = 


0.75) 


0.217 


± 


0.007 


0.40 


± 


0.04 


Model III (p = 


0.90) 


0.213 


± 


0.009 


0.39 


± 


0.04 



Table 1. Growth exponents of the isotropic on-lattice growth models. 



5. Summary 

In conclusion, we drew new considerations about the interface scaling of round clusters 
using large scale simulations of isotropic radial clusters, mainly the off-lattice Eden 
model. This approach reveals a subtle pitfall that can be present in the analysis of 
experiments with radial symmetry. Indeed, we show that the growth exponent depends 
on the strategy adopted to measure the interface width. For the particular case of 
the off-lattice Eden model, the expected value 1/3 is found when a fixed origin is 
used as reference. Otherwise, when the border CM is used as the origin, we found 
the growth exponent (3cm = 0.404 ± 0.013 very close to 2/5, in very good agreement 
with previous reports [23J. These differences arise from the border CM fluctuations 
increasing faster than those of the interface. We also show that the exponents (3q and 
Pcm are associated to the large and small wavelength fluctuations, respectively. These 
analyzes were corroborated by three distinct on-lattice models for which the lattice 
induced anisotropy is absent. 

It is important to emphasize that the essential features presented in this work were 
obtained from a very simple approach (inclusion of the CM analysis) not previously 
considered, while the models were used only as a support to these analysis. Moreover, 
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r 



Figure 6. Interface width evolution for the model III with p = 0.90. The averages 
were done over 140 independent samples, f and w are given in cell diameter unities. 

the difficulty for observing these features in the radial clusters previously studied, 
particularly the Eden model, lies on the lattice-induced cluster misshape that masks 
these fluctuations. Finally, it is also important to stress that the growth exponent 
has been used to draw conclusions about the underlying dynamics of systems with 
unquestionable scientific relevance (e.g. Refs. [9l [101 [11]) without caution with the 
pitfalls presented in this work. 
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